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Abstract — This paper attempts to answer a question: for a given 
traversal area, how to quantify the geometric impact of anchor 
placement on localization performance. We present a theoretical 
framework for quantifying the anchor placement impact. An 
experimental study, as well as the field test using a UWB ranging 
technology, is presented. These experimental results validate the 
theoretical analysis. As a byproduct, we propose a two-phase 
localization method (TPLM) and show that TPLM outperforms 
the least-square method in localization accuracy by a huge margin. 
TPLM performs much faster than the gradient descent method 
and slightly better than the gradient descent method in localization 
accuracy. Our field test suggests that TPLM is more robust against 
noise than the least-square and gradient descent methods. 

1. Introduction 

Accurate localization is essential for a wide range of appli- 
cations such as mobile ad hoc networking, cognitive radio and 
robotics. The localization problem has many variants that reflect 
the diversity of operational environments. In open areas, GPS 
has been considered to be the localization choice. Despite its 
ubiquity and popularity, GPS has its Achilles' heel that limits 
its application scope under certain circumstances: GPS typically 
does not work in indoor environments, and the power consump- 
tion of GPS receiver is a major hindrance that precludes GPS 
applications in resource-constrained sensor networks. 

To circumvent the limitations of GPS, acoustic/radio-strength 
and Ultra-wideband (UWB) ranging technologies have been 
proposed [10|, l22l . ifTTI . Measurement technologies include 
(a) RSS-based (received-signal-strength), (b) TOA-based (time- 
of-arrival), and (c) AoA-based (angle-of-arrival). 

RSS-based (received-signal-strength) technologies are most 
popular mainly due to the ubiquity of WiFi and cellular 
networks. The basic idea is to translate RSS into distance 
estimates. The performance of RSS-based measurements are 
environmentally dependent, due to shadowing and multipath 
effects (22|, E2, 0. TOA-based technologies such as UWB 
and GPS measures the propagation-induced time delay between 
a transmitter and a receiver. The hallmark of TOA-based 
technologies is the receiver's ability to accurately deduce the 
arrival time of the line-of-sight (LOS) signals. TOA-based 
approaches excel RSS-based ones in both the ranging accuracy 
and reliability. In particular, UWB-based ranging technologies 
appear to be promising for indoor positioning as UWB signal 
can penetrate most building materials ifTOl . 11221 . IfTTI . AoA- 
based technologies measure angles of the target node perceived 
by anchors by means of an antenna array using either RRS or 
TOA measurement [22], IfTTI . ||2- Thus AoA-based approaches 
are viewed as a variant of RRS or TOA technologies. 



Rapid advances in IC fabrication and RF technologies make 
possible the deployment of large scale power-efficient sensor 
networks. The network localization problem arises from such 
needs [26], [6|. The goal is to locate all nodes in the network in 
which only a small number of anchors know their precise posi- 
tions initially. A sensor node (with initially unknown position) 
measures its distances to three anchors, and then determines 
its position. Once the position of a node is determined, then 
the node becomes a new anchor. The network localization has 
two major challenges: 1) cascading error accumulation and 
2) insufficient number of initial anchors and initially skewed 
anchor distribution. The first challenge is addressed by utiliz- 
ing optimization techniques to smooth out error distribution. 
The second challenge is addressed by using multihop ranging 
estimation techniques [28], [27|. 

Centralized optimization techniques for solving network lo- 
calization include semidefinite programming by So and Ye 
lf24ll and second order cone programming (SOCP) relaxation 
by Tseng [26]. Local optimization techniques are realistic in 
practical settings. However, they could induce flip ambiguity 
that deviates significantly from the ground truth (6). To deal 
with flip ambiguity, Eren et. al. ||9] ingeniously applied graph 
rigidity theory to establish the unique localizability condition. 
They showed that a network can be uniquely localizable iff 
its grounded graph is globally rigid. The robust quadrilaterals 
algorithm by Moore et al. 12T1 achieves the network localiz- 
ability by gluing locally obtained quadrilaterals, thus effectively 
reducing the likelihood of flip ambiguities. Kannan et al. |fl2l . 
lfT3ll formulate flip ambiguity problem to facilitate robust sensor 
network localization. 

Lederer et al. |fT5l . [16] and Priyantha et al. If23l revolve 
around anchor-free localization problem in which none of the 
nodes know their positions. The goal is to construct a global 
network layout by using network connectivity. They develop an 
algorithm for constructing a globally rigid Delaunay complex 
for localization of a large sensor network with complex shape. 
Bruck, Gao and Jiang (6) establish the condition of network 
localization using the local angle information. They show that 
embedding a unit disk graph is NP-hard even when the angles 
between adjacent edges are given. 

Patwari et al. If22l used the Cramer-Rao bound (CRB) to 
establish performance bounds for localizing stationary nodes in 
sensor networks under different path-loss exponents. Dulman 
et al. El propose an iteration algorithm for the placement 
of three anchors for a given set of stationary nodes: in each 
iteration, the new position of one chosen anchor is computed 
using the noise-resilience metric. Bishop et al. 12, ID, 0, 12 



study the geometric impact of anchor placement with respect 
to one stationary node. Bulusu et al. proposed adaptive anchor 
placement methods [7|. Based on actual localization error at 
different places in the region, their algorithms can empirically 
determine good places to deploy additional anchors. 

Our work primarily focuses on the geometric impact quan- 
tification of an anchor placement over a traversal area. As a 
byproduct, it also forms the basis for optimal anchor selection 
for mitigating the impact of measurement noise. To our best 
knowledge, only a few papers 0, Q, El, E2, 0, 0, 
mentioned about the geometric effect of anchor placement but 
under a different context from this paper. The idea of quantify- 
ing the geometric impact of anchor placement on localization 
accuracy over a traversal area has not been considered before. 

The rest of this paper is structured as follows. Section 2 
shows examples to illustrate the anchor placement impact 
on localization accuracy. Section 3 proposes a method for 
quantifying the anchor placement impact. Section 4 presents a 
two-phase localization method Section 5 conducts a simulation 
study to validate the theoretical results. Section 6 presents 
the field study using the UWB ranging technology. Section 7 
concludes this paper. 

2. Effect of Anchor Placement 

Throughout this paper we use the term anchor to denote 
a node with known position. The positions of anchors are 
stationary and available to each mobile node (MN hereafter) 
in a traversal area l22l . The goal is to establish the position of 
a MN through ranging measurements to available anchors. The 
problem is formulated as follows: Let pi = (xi,yi), (1 < i < 
m) denote the known position of the i th anchor, and p the actual 
position of the MN of interest. The distance between p and pi is 
thus expressed as di — d(p : pi) — y/(x — Xi) 2 + (y — yt) 2 . In 
practical terms, the obtained distance measurement is affected 
by measurement noises. Thus the obtained distance is as 
di = di+e-i where di the true distance, and is widely assumed 
to be a Gaussian noise with zero mean and variance of |30|, 
), El, OH, 0, (22), UU. Define a function 
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where m is the number of accessible anchors, (xi,yi) is the 
position of the i th anchor, and di noisy ranging between the i th 
anchor and the MN. The localization problem ||22), 0, fl4) is 
formulated as 



min f(x,y),$t is the real number set. (1) 

(a;,y)eSR 2 

([T]) presents a nonlinear optimization problem that can be solved 
by a number of algorithms. This paper uses the gradient descent 
method (GDM hereafter) for its simplicity [1]. GDM is based 
on an intuitive idea that if f(x^°\ y(°>) is differentiable in the 
vicinty of (x^-°\y^), then f((x^-°\y^) decreases fastest in 
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where r\ refers to the iteration step size, i denotes the ith 
iterations. One starts with an initial value (x^-°\y^), then 
applies (j2j) iteratively to reach a local minimum (x,y), i.e., 

f(x,y)<f(x {l \y (l) ), l<i<n. 

There are three issues associated with GDM: 1) initial value; 
2) iteration step size and 3) convergence rate. The convergence 
rate could be sensitive to the initial value as well as the iteration 
step size. While in general there is no theoretical guidance 
for selecting the initial value in practice, the value obtained 
by a linearized method is often chosen heuristically as the 
initial value. This nonlinear optimization problem in ([1} can 
be linearized by the least-square method lTT4l . |fl9l as follows: 
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A T is the transpose of A, and = \J x 2 + y 2 , 1 < i < m. The 
least-square method (LSM hereafter) uses noisy measurements 
to estimate the position of MN 




Fig. 1. Hilbert Trajectory 

To shed light on the anchor placement effect on localization 
accuracy, consider a Hilbert traversal trajectory with different 
anchor placement (AP) setups in Table [T] The Hilbert trajectory 
is formed by the piecewise connection of 8190 points in an 
100 x 100 region (see Fig((T|i). A Hilbert trajectory (HT here- 
afer) has a space-filling property. As a result, the localization 
performance obtained on a HT is a good approximation to that 
of the underlying area. 

We conduct an experiment as follows: start from the upper- 
left corner, the MN moves along the HT. At each point p (the 
position of MN), noisy distance measurements di = di + e; 
from p to the anchors are generated. Both LSM and GDM are 
then used to derive the estimated position. To avoid artifacts, 



for a given AP and noise level, the error statistics in Table [T] 
are obtained by traversing the HT 10 times. Each HT traversal 
involves the establishment of 8190 positions using both LSM 
and GDM (GDM uses the estimated position derived from LSM 
as its iteration initial value), with a step size of 0.00001. The 
termination condition of GDM is set as |[p( l+1 )— p^||2 < 0.001 
and the maximum number of iterations is set as 100. 

TABLE 1 

Anchor placement & localization accuracy 
Least-Square Method 



positions of anchors 


a 


ave 


std 


time 


(0,100), (0,0), (100,0) 


0.3 


0.42 


0.26 


2.52 


1.0 


1.40 


0.85 


2.54 


(0,100), (7, 50), (3, 40) 


0.3 


4.01 


3.65 


2.51 


1.0 


13.45 


12.15 


2.57 


Gradient Descent Method 


positions of anchors 


(7 


ave 


std 


time 


(0,100), (0,0), (100,0) 


0.3 


0.37 


0.21 


63 


1.0 


1.22 


0.69 


81 


(0,100), (7, 50), (3, 40) 


0.3 


0.76 


0.92 


275 


1.0 


2.94 


7.41 


293 



Gaussian noise JV(0, a 1 ) 



The ave and std fields in Table [T] denote average local- 
ization error and standard deviation, and the time field the 
elapsed time per HT traversal (seconds). Table [T] shows that 
AP (0, 100), (0,0), (100,0) yields a much better localization 
accuracy than AP (0, 100), (7,50), (3,40) for both LSM and 
GDM. This indicates, a fortiori, that an anchor placement (AP) 
has a significant bearing on localization accuracy over an area. 
This observation raises a question: can we quantify the impact 
of an anchor placement (AP) with respect to an area? The aim 
of this paper is to answer this question. 

3. Anchor Placement Effect Quantification 

The section focuses on the effect quantification of anchor 
placement. Our main device is based on the notion of the 
geometric dilution of precision (GDOP) ED, fl/7), (22), (8). 
We begin with the anchor-pair GDOP function as follows: 

Theorem 1: Let pi and Pj be the positions of anchor pair 
i and j, and di and dj be the actual distances from MN at 
p = (x,y) to the anchors then the geometric dilution of 
precision at p with respect to Pi,Pj, denoted by the anchor-pair 
GDP function g2(Pi,Pj)(p), is 
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where \\p~ Pj\\ is the distance between pi — (xi,yi) and pj 
Proof. Let h(pi,pj)(p) be a matrix defined as 
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where (sin(a), cos(a)) and (sin(/3), cos(/3)) denote the direc- 
tion cosines from p to anchors at pi and pj. The anchor-pair 
GDOP function <?2(Pi,Pj)(p) is 



92(Pi,Pj)(p) = \Jtr({h T h)- 
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where T/tr denotes the transpose/trace of a matrix, and A^ 1 
refers to the inversion of A. A simple manipulation obtains 

det(h T h) = sin 2 (/3 — a) and 
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into 

□ 



<|8J> yields 

Theorem"^] asserts that /3 — a degree separation represents 
the effect of anchor pair at Pi,Pj on the localization accu 
racy at p. It implies that j3 — a = tt/2 degree separation 
(32 (Pi j Pj ) (p) = v2) are best in localization accuracy, whereas 
(3 — a — degree separation (colinear) (g2(Pi,Pj)(p) = 00) are 
worst. The smaller the <72(Pi>Pj)(p)> the better the accuracy of 
localization. This is a well-known result also obtained in [22 1, 

a, 0, a m. 

We now extend the anchor-pair GDOP function into the 
multi-anchor GDOP function. Let pi , • • • , p m be a set of 
positions of anchors accessible to MN at p. The multi-anchor 
GDOP function g m (pi, ■ ■ • ,Pm)(p) is linked to the anchor-pair 
GDOP function as follows: 



t(Pl) ' ' ' ,Pm)(p) 



min g 2 (Pt,Pj)(p) 

1< z ,j<rn 



(9) 



We call an anchor pair (pi,Pj) the optimally selected 
anchor pair (OSAP) if 172 (Pi , Pj ) (p) has a minimum value 
among all anchor pairs from from p%, ■ ■ ■ ,p m . |9]i shows that 
3m (pi, • • • ,Pm)(p) is determined by the OSAP, which in turn 
varies with p and accessible anchors. 

Figs(|2|i-(j3|l visualize <?3(---)(p) function under the two 
anchor placements over the 100 x 100 region. A three- 
dimensional graph in Figs(|2]i-(|3]l, which is called the least 
vulnerability tomography (LVT), geometrizes the anchor place- 
ment effect: the LVT elevation has an implication: when 
in a trough area, the noise has less impact on localization 
accuracy than when in a peak area. FigQ shows that the 
LVT of AP (100,0), (0,0), (0,100) has terrain waves with 
a elevation variation from v2 to 2. In contrast, the LVT 
of AP (0,100), (7,50), (3,40) shown in Fig([3jl is relatively 
flat for the most part of the region but has a dramatic 
elevation variation from 2 to 16 in the vicinity of (0, 0) 
and (0, 100). Overall, the terrain waves in FigQ has a 
much higher elevation than that in Fig(j2j. This offers an 
explanation: why AP (0, 100), (0, 0), (100, 0) outperforms AP 
(0, 100), (7,50), (3,40) as shown in Table [l] 

The LVTs in FigsQ-(|3) induce us to extend the GDOP 
function from a point into an area. Let g m (pi,--- ,p m )(ft) 
denote the average elevation of g m (pi, ■■■ ,Pm)(p) over ft. The 




Fig. 2. LVT of anchor placement (0, 100), (0, 0), (100, 0) 




Fig. 3. LVT of anchor placement (0, 100), (7, 50), (3, 40) 

relation between g m (p u ■ ■ ■ ,p m )(fi) and g m {pi, ■ ■ ■ ,p m )(p) 
(g m (pi,--- ,Pm)(x,y)) becomes 

1 



9m(P 
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mi 



,Pm)(x,y)dxdy 
(10) 

Similarly, the notion of g m (pi, ■ ■ • ,Pm)(p) can be extended 
from a point p to a trajectory T as 

1 



9m(PU ■ ■ ■ )Pm)(T) = •TpT J g m (pi, ■ ■ ■ ,p m )(x,y)dl (11) 

The discrete form of ( fTT| becomes 
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(10 1 provides a means for quantifying the impact of an AP 
over an area. In practice, due to the arbitrariness of the area 
boundary and of an anchor placement, it would be impossible 
to derive a closed-form expression. In this paper, we use 
Trapezium rule method to compute ( fTO] ). It is done by first 
splitting the area into 10, 000 non-overlapping sub-areas, then 
applying the Trapezium rule on each of these sub-areas. 

TABLE 2 

Anchor Placement Impact Quantification 



Anchor Placement 


pi = (0, 100), pa = 


(0,0), p 3 = (100,0), p 4 = (7,50), p 5 = (3,40) 


Anchor Placement Impact over ft 


93(Pl,P2.P3)(n) 


93(Pl,P4,P5)(fi) 


1.501 


2.615 


Anchor Placement Impact over F 


93((Pi,Pa,P3)(r) 


93(Pl,P4,Ps)(r) 


1.499 


2.622 


Q: 100 X 100 areaj: HT 



Table [2] provides a numerical calculation of AP impact over 
Q and over T. The AP impact calculated in Table [2] implies 



that the localization accuracy over the AP\ is better than that 
under AP 2 as g 3 (p 1 ,p 2 ,p 3 ){n) < g 3 (p 1 ,p 4 ,P5){ty, which 
agree with the results obtained in noisy environments in Table [Tj 
Comparing g 3 (p 1 ,p 2 ,p 3 )(ty and g 3 (pi,P2,P3)(r) in Tableg 
shows that the difference between them is very small. However, 
computation of g 3 {pi,P2,P3)(ty takes about 23 seconds, as 
opposed to 0.2 seconds in computing g 3 {pi,p2,pz)iT)- 

4. Two-Phase Localization Algorithm 

It is clear from Fig(|2| that g m (pi, • • • ,Pm)(p) is determined 
by an OSAP and the OSAP varies by area. FigQ shows the 
OSAP areas under AP (0, 100), (0, 0), (100, 0). Observe the 
noise-free scenario in the left-hand side graph in FigQ where 
each colored area represents an anchor pair chosen over its two 
alternatives based on d9j. The white colored areas represent the 
regions in which the anchor pair at (0, 100), (0,0) is chosen. 
The blue colored areas represent the regions where the anchor 
pair at (0, 0), (100, 0) has its geometric advantage over its 
alternative anchor pairs. The red colored areas represent the 
regions where the anchor pair at (100, 0), (0, 100) is expected 
to produce the most accurate localization. The right-hand side 
graph in FigQ plots the OSAP areas under the noise level 
of a = 1.0. The presence of noise makes the borderlines on 
different colored areas rough and unsmooth. This is because 
that borderlines correspond to the isolines by two anchor pairs 
and hence are highly sensitive to noise. 





Fig. 4. a) noise free b) noise level a(1.0): white area by anchor pair 
at (0, 100),(0, 0), black area by anchor pair at (0, 0), (100, 0), red area 
by anchor pair at (0, 100),(100, 0) 

This observation of FigQ leads to a new two-phase localiza- 
tion algorithm (TPLM) that differs significantly from LSM and 
GDM. TMPL differs from LSM and GDM. It proceeds in two 
phases: In the first phase called optimal anchor pair selection, 
an optimal anchor pair at (pi,Pj) is identified based on (|9ji, 
then the position(s) can be estimated by solving 



Pi 



Pj\ 



(13) 



Through a lengthy manuplation, two possible solutions to ( 13 i 
are obtained as follows 
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In the second phase called disambiguation, the reference point 
obtained from LSM is used to single out one from two possible 



positions (14i. It is done by choosing one position that has a 
shorter distance to the reference point. 

Algorithm 1 Two-phase Localization Method 

Input: pi = (xi,yi), 1 < i < m: i th anchor's position 
Input: df. distance measurement from j th anchor to MN p 
Input: p = (x,y): a reference point obtained via LSM 



9min <~ oo,p+ <- (0,0), p ++ <- (0,0) 
for i = 1 to m — 1 do 
for j = i + 1 to m do 
if g2(puPj)(p) < 9mm then 

g mi n <- g(pi,Pj)(p), (p+,p++) <- (puPj) 
end if 
end for 
end for 

d = \p+ -P++\,d+ = \p-p+\,d+4 
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The code between lines |2|8] examines each anchor pair in 
turn and chooses an optimal anchor pair, which corresponds 
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constitute the basic 



to (|9jl. The complexity is (™). Lines 
block of the optimal anchor selection phase. The code between 
lines [T3p6] is used to disambiguate the two possible solutions 
using the reference position p obtained by LSM. 

TABLE 3 

Two-Phase Localization Method 



positions of anchors 


cr 


ave 


std 


time 


(0, 100), (0,0), (100,0) 


0.3 


0.40 


0.22 


3.17 


1.0 


1.32 


0.73 


3.17 


(0, 100), (7,50), (3,40) 


0.3 


0.70 


0.88 


3.44 


1.0 


2.76 


5.67 


3.15 


Gaussian noise JV(0, a 1 ) 



We present the results of TPLM in Table [3] using the same 
anchor placement setups in Table [T] In all cases, TPLM gives a 
significant error reduction over LSM with a minor degradation 
as TPLM uses the reference point obtained from LSM. It is 
fairly clear that TPLM is an order of magnitude faster than 
GDM, and performs slightly better than GDM in terms of 
localization accuracy. 



5. Simulation Study 

The aim of this section is threefold: first, to conduct the 
performance comparison between TPLM, LSM, and GDM 
in noise environments; second, to study the impact of noise 
models; and third, to investigate the random anchor placement 
impact. 

A. Visualizing Noise-induced Distortion 

To visualize the difference among LSM, GDM, and TPLM, 
we plot the restored HT by LSM, GDM, and TPLM under 
Gaussian noise model in Higs([5j-(j7j using AP4 in Table [4] 
Notice that anchor positions are plotted as solid black circles. 

TABLE 4 
Anchor Placement Setup 



Anchor Position 


Pi = 


(0,100),p 2 = (7,50),p 3 = (0,0) 


Pi = (3,40),p 5 = (100,0),p 6 = (1,98) 


Anchor Placement Setup 


APi 


(0, 100), (0,0), (100,0) 


AP 2 


(0,100), (7, 50), (3, 40) 


AP 3 


(0, 100), (0,0), (100,0), (1,98) 


AP A 


(0, 100), (7, 50), (3,40), (1,98) 





(a) (b) 

Fig. 5. restored HT by LSM: a) noise (0.02); b) noise (0.2) 





(a) <b) 

Fig. 6. restored HT by GDM: a) noise (0.02); b) noise (0.2) 

Visual inspection reveals the apparent perceived difference 
among the restored HTs by LSM, GDM and TPLM: when 
the noise level a is 0.2, the restored HT by LSM becomes 
completely unrecognizable in Fig{5])(b). While in Fig([6])(b) the 
upper right portion of the restored HT by GDM to some degree 
preserves the hallmark of the HT, but the most part of the 
recovered HT is severely distorted and barely recognizable. 
The restored curve by TPLM contrasts sharply with those 



Fig. 7. restored HT by TPLM: a) noise (0.02); b) noise (0.2) 



by LSM and GDM in its preservation of fine details of the 
HT for the most part, while having minor distortion in the 
lower- and upper-left corner areas in Fig|7|(b). Such spatially 
uneven localization performance can be explained by examining 
FigQ, which shows that the lower- and upper-left areas exactly 
correspond to the LVT peak areas. The perceived differences 
between GDM and TPLM in Figs|5])-(|7]i can be quantified as 
follows: GDM has the average error of 0.69 per HT traversal, 
while TPLM produces the average error of 0.465. In addition, 
GDM takes about 218.48 seconds per HT traversal, in contrast 
to 2.94 seconds taken by TPLM. 

B. Gaussian Noise vs. Non-Gaussian Noise 

In this subsection, we will compare the performance of GDM 
and TPLM under Gaussian and uniform noise models, there are 
some scenarios where ranging noise may follow uniform model 
0, EDI . In the experimental study, the performance between 
GDM and TPLM is compared under a same noise level with the 
different noise models. For Gaussian noise model A/"(0, a 2 ), a 
refers to the noise level in Fig([8]l, and for a uniform noise model 
U(—a, a), the noise level in Fig(|9jl is expressed as a 2 /3. The 
data presented reflects the average error of GDM and TPLM 
over 10 HT traversals under the anchor placements in Table [2] 




gaussian noise: noise level 

Fig. 8. Accuracy: (top) GDM; (bottom) TPLM (Gaussian noise) 

The graphs in Figs([8j-([9j represent the localization error 
curves of GDM and of TPLM, respectively. Under both the 




uniform noise: noise level 



Fig. 9. Accuracy: (top) GDM; (bottom) TPLM (uniform noise) 



noise models, TPLM outperforms GDM by huge margins under 
AP 2 and AP 4 . Both GDM and TPLM perform indistinguish- 
ably under AP\ and AP% . While the impact difference between 
APi and AP% is barely noticed as their performance curves 
are overlapped, the performance curves between AP\ and APi 
and between AP 3 and AP4 are clearly separated, thus the 
impact difference AP\ and AP2 and between AP3 and AP4 
are evident. 
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Fig. 10. Elapsed time: (top) GDM; (bottom) TPLM (normal noise) 

The results show that the Gaussian noise has more impact on 
GDM than the uniform one. By comparison, the performance of 



TPLM is insensitive to noise models. Fig( 10 1 plots the elapsed 



time by GDM and TPLM per HT traversal, showing that TPLM 
is orders of magnitude faster than GDM. 

C. Random Anchor Placement 

We study the impact of random anchor placement on the 
performance of LSM, GDM and TLPM. To achieve this, the 
number of anchors are randomly placed in 100 x 100 and 
50 x 100 regions. For each randomly generated anchor place- 
ment (RGAP) in a region, the g m (pi,- ■■ ,p m )(T) function is 
calculated, and the error statistics of LSM, GDM, and TPLM 
under Gaussian noise level of 0.3 are gathered and compared. 



Fig( 11 i(a)-(b) show the actual g m (pi, ■ ■ ■ ,p m )(T) (average 
LVT elevation over T)) distribution by 100 RGAPs. In the 
top graph, 3 anchors are randomly placed (r.u.p.) in the entire 
traversal area while in the bottom graph 3 anchors are r.u.p. in 
the upper half of the traversal area. Fig([TT]i(a) clearly exhibit 
positive skewness. This implies that a RGAP over the entire 
traversal area in general yields good localization accuracy. 
Recall that the smaller g m (pi, ■ ■ ■ ,p m )(T) is, the better local- 
ization is. The g m {pir'' >Pm)(r) distribution in Figdli(b) 
is more dispersed than that in the Fig([TT}(a), meaning that a 
RGAP over the half traversal area in general underperforms a 
RGAP over the entire traversal area. 




Fig. 11. (a) g 3 (...)(T) distribution by 100 RGAPs with 3 anchors in 
100 x 100; (b) g 3 (...)(T) distribution by 100 RGAPs with 3 anchors 
in 50 x 100 



TABLE 5 
Random Anchor placement 

Random Anchor Placement in The Traversal Area 



Anchors 


3m (Pi, •■Pm)(r) 


TPLM 


LSM 


GDM 


ave 


ave 


ave 


3 


2.05 


0.57 


1.52 


0.87 


4 


1.69 


0.44 


0.79 


0.49 


5 


1.53 


0.40 


0.59 


0.38 


6 


1.48 


0.39 


0.51 


0.33 


Random Anchor Placement in the Half Traversal Area 


Anchors 


9m(pi, ..Pm)(r) 


TPLM 


LSM 


GDM 


ave 


ave 


ave 


3 


2.72 


0.97 


2.82 


1.54 


4 


1.93 


0.52 


1.37 


0.69 


5 


1.77 


0.46 


0.81 


0.46 


6 


1.69 


0.44 


0.72 


0.40 



Fig( 12 1 presents the performance curves of LSM, GDM, and 
TPLM under RGAPs, where the x axis refers to the value of 
gm(pi, •■Pm)(r) and the y axis the error deviation and average 
error. As is clear from Fig ([12}, TPLM outperforms both LSM 
and GDM in localization accuracy and reliability: both the error 
deviation and average error curves of TPLM are consistently 
lower than those of LSM and GDM. 

Table [5] tabulates the results of RGAPs with different num- 
ber of anchors, where g m (pi, ..,p m )(T) denotes the average 
3m(pi, •■Pm)(r) over 100 RGAPs with a fixed number of 
anchors. It shows that the localization accuracy of TPLM over 
GDM is diminished as the number of anchors increases. It 
becomes evident that under the noise level of 0.3, the actual 




Fig. 12. (a) error deviation; (b) average error (100 RGAPs with 3 
anchors over the 100 x 100 region, a — 0.3) 



performance of LSM, GDM, and TPLM deteriorates as the 
value of g m (pi, ..,Pm)(r) increases. This solidifies the fact 
that g m (pi, •■Pm)(r) is an effective discriminator for the anchor 
placement impact on localization performance. 

6. Field Experimental Study 

This section focuses on the field test of a DARPA-sponsored 
research project, using the UWB-based ranging technology 
from Multispectrum Solutions (MSSI) iflOl . and Trimble dif- 
ferential GPS (DGPS). Our field testbed area was a 100 x 100 
square meters. It consisted of an outdoor space largely oc- 
cupied by surface parking lots and an indoor space inside a 
warehouse. This testbed area was further divided into 10, 000 
non-overlapping lxl grids. Each grid represents the finest 
positioning resolution to evaluate RF signal variation. 

TABLE 6 

Field testbed area and anchor placement 





V'o 




Vf 


a 


-74.476069 


40.537808 


84719 


111045 


0.381583 


Anchor placement 


GPS position 


transformed position 


-74.475585 


40.538468 


p 4 = (65.345179,52.75145) 


-74.475287 


40.53856 


p 5 = (92.580022,52.83239) 


-74.475186 


40.538294 


P6 = (* 


59.52274, 22.232383) 



In the field test, DGPS devices were mainly used for outdoor 
positioning while MSSI UWB-based devices were used for 
indoor positioning. The experimental system was composed of 
four MNs equipped with both the Trimble DGPS and MSSI 
UWB-based ranging devices. Three UWB devices were used 
as anchors being placed at known fixed positions. Using the 
known positions of the UWB anchors and real-time distance 
measurements, each MN could establish the current position as 
well as that of other MNs locally, at a rate of approximately 2 
samples per second. Each MN could also establish the position 
via the DGPS unit at a rate of 1 sample per second when 
residing in the outdoor area. The field test area was divided 
into four slightly overlapping traversal quadrants denoted as 
(r ,ri,r 2 ,r 3 ). A laptop (MN) equipped with DGPS and 



UWB ranging devices was placed in a modified stroller as 



shown in Fig(13i, and one tester pushed the stroller, traveling 



along a predefined trajectory inside a quadrant (see Fig(14i) 
Each run lasted about 30 minutes of walk. 



Rain Cover. 



Basket w/GPS 
unit & battel i 




GPS Antenna 



UWB radio 



Laptop running 
system s/w 



Fig. 13. Stroller with a laptop attached with DGPS and UWB ranging 
devices 




Fig. 14. Four trajectories (Fo, Ti,T2, ^3), black circles: the positions 
of anchors, yellow rectangle: warehouse 

With repeated trial runs, we found that the UWB signal could 
barely penetrate one cement wall of the warehouse. To establish 
positions inside the warehouse, we placed one anchor close to 
the main entrance of the warehouse while placing two other 



anchors around the center of the field testbed area. Fig( 14 1 
showed the four trajectories (r , r 1; r 2 , T 3 ) in the field test, 
which are translated from the GPS trajectories. A part of Tq 
trajectory was inside the warehouse, thus the part of GPS of 
To was not available. Table|7] provides the calculated results 
°f 33(P4>P5>P6)(r) using the GPS trajectory data, indicating 
that r would yield the most accurate localization while Fx 
produced the least accurate localization. 

TABLE 7 

33(P4iP5;P6)(r) OF TRAJECTORIES IN FIELD TESTBED 



33(P4,P5,P6)(ro) 


S3(P4,P5,P6)(ri) 


1.667025 


1.967605 


33(P4,P5,P6)(r 2 ) 


93(P4,P5,P6)(r3) 


1.787370 


1.835048 
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Fig. 16. Restored Trajectories by (a) TPLM, (b) GDM, (c) LSM (black 
circles: positions of UWB anchors) 



In the field testbed, the trajectory of each MN was controlled 
by an individual tester during a 30-minute walk. For purpose 
of the primary experiment being conducted, the strollers were 
fitted with bicycle speedometers to allow the tester to control 
his speed, in order to produce sufficient reproducibility for the 
primary experiment, but not for our testing of localization. Due 
to inherent variability in each individual movement, an objec- 
tive assessment of localization accuracy without a ground truth 
reference is almost impossible. For a performance comparison, 
we used GPS trajectories as a reference for visual inspection 
of the restored trajectories by LSM, GDM, and TPLM. 



Three curves in Fig( 15 i(a)-(d) represent the field distance 
measurements between a tester and the three UWB anchor de- 
vices during a 30-minute walk in trajectories (To, Ti, T2, T3). it 
is fairly obvious that measurement noise in different trajectories 
vary widely: The distance measurement curves in trajectories 
(r ,r 3 ) are discontinuous and jumpy in Figs( 15 l(a)&(d), in 



contrast to the relatively smooth distance curves in trajectories 
(Ti,r2) in Fig([T5]l(b)&(c). Such a discontinuity in measure- 
ments occurred when testers were traveling in an area where 
UWB devices on stroller had no direct line of sight to UWB 
anchors, thus introducing additional noise in r ,r 3 . 



Figs( 16 i(a)-(c) show the restored trajectories by TPLM, 
GDM, and LSM using the UWB ranging technology in the 
field test. A visual inspection suggests that LSM was extremely 
prone to noise as its restored trajectories were appreciably 
distorted beyond recognition in some parts. As indicated in 
Fig(fl6|(b), GDM gave an obvious error reduction over LSM 
but at the expense of computational cost. The most visually 
perceived difference between LSM and GDM can be seen in 
the circled areas in Figs([T6|(b)-(c). By contrast, the difference 
between GDM and TPLM can be visualized in the circled areas 



Figs( 16 i(a)-(b) where TPLM produced a detail-preserved but 
slightly distorted contour of the trajectory. A further offline 
analysis showed that on average GDM takes 18.87ms per 
position establishment, while TPLM/LSM take 0.4/0. 35ms. 




Fig. 15. Distance measurements to three UWB anchors: (a) To, (b) fi, (c) T2, (d) T3 



7. Conclusion 

This paper studies the geometric effect of anchor placement 
on localization performance. The proposed approach allows the 
construction of the least vulnerability tomography (LVT) for 
comparing the geometric impact of different anchor placements. 
As a byproduct, we propose a two-phase localization method. 

To validate theoretical results, we conduct comprehensive 
simulation experiments based on randomly generated anchor 
placements and different noise models. The experimental results 
agree with our anchor placement impact analysis. In addition, 
we show that TPLM outperforms LSM by a huge margin in 
accuracy. TPLM performs much faster than GDM and slightly 
better than GDM in localization accuracy. The field study shows 
that TPLM is more robust against noise than LSM and GDM. 

In this paper, we adopt a widely used assumption that 
measurement noise is independent of distance. This assumption 
is somewhat unrealistic in some practical settings. In our future 
research, we will perform experiments using UWB ranging 
technology to study the noise-distance relation in indoor and 
outdoor environments. The anchor placement impact analysis 
will be refined to incorporate more realistic noise-distance 
models. 
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